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Motivation 

Penetration and absorption of solar radiation in the upper layers of the ocean is affected by 
the optical properties of the water column. The optical properties of the water column in turn 
are affected by its biochemical composition. The current NASA Ocean Biogeochemical 
Model (Gregg, 2008) is coupled with an OGCM and forced by prescribed atmospheric data. 
In anticipation of the integration of the biochemical model into the next generation GMAO 
fully coupled system, an intermediate solution has been implemented to estimate the 
penetration depth {l/Kd_PAR) of ocean radiation based on chlorophyll concentration, thus 
allowing for the feedback between the ocean biology component and the general circulation 
model. 

Data and experiments 

The GEOS-5 coupled ocean-atmosphere-sea ice model was used for this study. Its 
components are: the MOMS ocean model running at 1/2° horizontal resolution with 40 
vertical levels, the GEOS-5 atmospheric model running at 1.5°xl° horizontal resolution with 
72 vertical layers and the sea ice model LANL CICE. The chlorophyll is modeled as a 
tracer with sources/sinks coming from the assimilation of MODIS chlorophyll data. The 
univariate SAFE (Keppenne et all 2014) was used to processed the chlorophyll data. 

Two experiments were conducted. In the first, CHLO , climatological values of Kd_PAR 

were used. In the second, CHLO , observed daily chlorophyll concentrations were 
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assimilated and Kd_PAR was derived according to Morel et al (2007). No other data was 
assimilated to isolate the effects of the time-evolving chlorophyll field. 

Figure 1 illustrates the amount of data available during a single assimilation cycle, while 
figure 2 shows the observation count for an entire year. The high latitudes and areas of 
persistent cloud cover are deprived of observations. 
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Figure 1: A5-(days composite of chlorophyll 
concentration anomaly with respect to the 1998-2006 
climatology. 
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Figure 2: Chlorophyll observations count for 1 year; the 
observations are binned to the ocean model grid. 
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Evaluation 

The daily MODIS Kd_PAR product was used to validate the skill of the penetration depth 
estimation and the MERRA-OCEAN re-analysis was used as a benchmark to study the 
sensitivity of the upper ocean heat content and vertical temperature distribution to the 
chlorophyll input. The metrics used to compare the results are the absolute error and the 
normalized root mean square error, expressed in terms of percentage. In the formula below, 
X stands for CHLO , Y for CHLO and Z for MERRA-OCEAN. 

dim assim 


Root Mean Square Error = X :Z = 
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Normalized Root Mean Square Error = 


^RMSE ^ ^ ^RMSE ^ ^ 
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Penetration depth validation 
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Figure 3: The difference of the absolute errors in 
the penetration depth with respect to MODIS 
Kd_PAR. 


Figure 3 shows a sample 5-day snapshot 
of the penetration depth from observations 
compared to simulations. The data 
coverage is spotty, so to obtain a more 
robust measure of this difference, a spatial 
integral of the error was computed and is 
shown in figure 4. Here the global integral 
is uniformly smaller for CHLO by 

about 1 m, similar to the integral over 
tropical Pacific Ocean (30°S-30°N). The 
error is generally larger for both 
experiments and exhibits stronger 
seasonal variability in the tropical Atlantic 
Ocean. 



Figure 4: Time series of the integrated penetration depth error, 
CHLO minus OBS is shown by dashed lines, CHLO minus 

OBS by solid lines. 


A similar seasonal 
enhancement in error 
amplitude variability in 
CHLO is observed in the 

assim 

Indian Ocean, which reflects 
the interannual variability in 
chlorophyll concentration. The 
overall reduction in 
penetration depth error may in 
part be attributed to the 
interdecadal variability of 
chlorophyll concentration, 
since the reference 
climatology period does not 
overlap with the experiment 
dates. 


Upper ocean heat content 

Figure 5 shows the map of A of the 

mean temperature of the upper 50 meters. 
The pink color indicates the smaller error 
in CHLO over CHLO . The most 

assim dim 

uniform improvement is observed in the 
Pacific basin. In the Atlantic, CHLO 

assim 

shows better performance only in the 
areas off the equator . The error 
difference in the Indian Ocean is generally 
smaller, being slightly in favor of CHLO 

in the east and CHLO in the west. 

dim 
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Figure 5. A for the mean upper 50 m temperature. 

^ NRMSE ^ 


Temperature vertical profiles 
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Figure 6. Map of stationary moorings in the 
Pacific and Atlantic oceans. The red squares 
mark the locations presented in this study. 


The maps in figure 6 highlight the 
regions most responsive to chlorophyll 
variability. To look at the effects of the 
chlorophyll assimilation on the vertical 
structure of the water column, we 
selected several locations coinciding 
with the stationary moorings of the 
TAO/TRITON and PIRATA arrays. 

These locations are along the equator in 
the Pacific Ocean and in the eastern 
Pacific Ocean cold tongue, as well as in 
the tropical Atlantic Ocean off the east 
coast of South America. 

Figure 7 displays the difference in 
absolute error with respect to the 
MERRA-OCEAN reanalysis at the 
mooring locations along the equator in 
the Pacific Ocean. Blue and purple solid 
lines show the absolute error for 
CHLO and CHLO , respectively. 

Light pink and blue shading show A 

^ ^ ^ NRMSE 


Pink shading corresponds to improvement in matching the reanalysis temperature, i.e. 
correcting the model bias, in CHLO over CHLO . The positive impact of assimilating 

chlorophyll concentration propagates downward to 60 meters in the east and 200 meters 
in the west. This roughly corresponds to the thermocline depth. The relative gain in skill 
is between 2% to 5%. Looking along the 95°W meridian (figure 8), one can see the 
deeper impact of the chlorophyll assimilation off the equator than at (95°W, 0°). 


As one can see from map of the upper ocean heat content A (figure 5), there is no 

significant improvement in the equatorial Atlantic Ocean in the CHLO experiment 

over the CHLO run. The positive impact of chlorophyll assimilation on <T(50m)> is 

found between 30°S and 10°S. Thus, we selected moorings located in this region to 
analyze the vertical temperature profiles (these are marked by red squares in the bottom 
plot on figure 6). The results are mixed, with the 8°S location showing a small 
improvement in skill (1%) in the upper 40 m and then a negative impact below (up to 
6%). The 14°S and 19°S locations both show the positive impact of the chlorophyll 
assimilation (5-12%) extending down to 120-150 meters. 
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Figure 7. 


Figures 7-9: The purple line 
shows CHLO , , the blue 

line shows CHLO the 

assim 

absolute temperature error 
computecj with respect to 
the MERRA-OCEAN 
reanalysis which 
assimilates all available 
in-situ observations, 
including TAO/TRITON and 
PIRATA moorings. The 
corresponding scale is on 
the bottom. The scale on 
the top shows A in 

percent with areas of 
positive values shaded in 
pink and negative values 
shaded in light blue. 
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Figure 8. 


Conclusion 

In the experiment with daily chlorophyll data 
assimilation, the penetration depth was estimated 
more accurately, especially in the tropics. As a result, 
the temperature bias of the model was reduced. A 
notably robust albeit small (2-5 percent) improvement 
was found across the equatorial Pacific ocean, which 
is a critical region for seasonal to interannual 
prediction. 

This study highlighted problems with simple 
chlorophyll assimilation, such as poor definition of the 
chlorophyll concentration in poorly observed regions. 
The implementation of the full biogeochemical model 
may partly alleviate these problems. 

The current limited study will be compared to the 
assimilation of chlorophyll in a full ocean 
biogeochemical model. 
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Figure 9. 





